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We have investigated the processes N(7r, 7r)N and N(7r, 77)N close to eta threshold using a simple, 
nonrelativistic Lee model which has the advantage of being analytically solvable. It is then possible 
00 ' to study the Riemann sheets of the S-matrix and the behavior of its resonance poles especially 

On . close to threshold. A theoretical simulation of the experimental cusp effect at eta threshold leads 

to a characteristic distribution of poles on the Riemann sheets. We find a pole located in the 4 
Riemann sheet that up to now has not been discussed. It belongs to the cusp peak at eta threshold. In 
, addition we obtain the surprising result using the Lee model that the resonance 5*11(1535) does not 

D ' play a large role. The main features of the experimental data can be reproduced without explicitly 

pL^ I introducing this resonance. Furthermore, we have also studied the reactions N(7, 7r)N and N(7, r;)N 

. and find reasonable agreement between the data and both models with and without the 511(1535) 

(N 



X 



resonance. 
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O . I. INTRODUCTION 

(N 

' The advent of new electron accelerators and intense photon sources substantially improved the data basis of electron 
Q.^ , and photoproduction of eta mesons. High quality data for angular distributions and total cross sections for photon 
■~»>^ ■ energies between threshold and 790 MeV may be considered a qualitative breakthrough in the experimental field. The 
[ progress in experimental proficiency demands a better theoretical understanding of these processes. 
I . The eta meson is believed to be mainly produced by the resonance iS'ii(1535) which is located close to eta threshold 
C_) (1486 MeV). Considerable work has been done in the case of eta photoproduction [Q, ||^, with the 511(1535) 
^ resonance as the main production channel. 

■ The eta threshold also infiuences scattering cross sections other than eta production. This is especially true in 
^ , elastic pion scattering where a strong cusp effect is observed due to unitarity. The presence of this cusp makes it 
difficult to get information on the resonance S'ii(1535), because the cusp and the resonance have similar signatures 
when analyzing experimental data. Even the existence of the resonance 6*11 (1535) seems to be questionable Q], Q|. In 
^ . a recent coupled channel calculation for eta and kaon photoproduction, Kaiser et al. Q have shown that the 511(1535) 
can as well be explained as a quasi-bound state of kaon/S-hyperon. 

It is generally accepted that the 511(1535) couples strongly to the eta while the neighboring resonance 5ii(1650) 
for all practical purposes does not decay into eta mesons. It is not clear why these two neighboring resonances, having 
the same quantum numbers, behave so differently. 

We investigated theoretically the combination of a resonance and the threshold effect using a Lee model because of 
its simplicity and clearness. It is fully analytically solvable, unitary and can in principle be extended to an arbitrary 
number of different mesons and resonances. The absence of antiparticles makes the model nonrelativistic. However, a 
covariant formalism of resonance excitation also leads to problems, e.g. anti-resonances (qqq) that are very unlikely 
in terms of quark models. Analyticity is mainly due to an inherent Tamm-Dancoff approximation in the model which 
limits the number of mesons present at any instant. 

Our model describes the interaction of the nucleon, the two resonances 511(1535) and 5ii(1650), the pion and the 
eta meson. Making use of perturbation theory we also included photoproduction of pions and eta mesons. From our 
calculation we find that two-pion production can not entirely be neglected. It was accounted for in a simple manner 



as discussed in the sect. [H. Coupling constants used in our model were determined by fitting to experimental data 
from pion scattering experiments and eta (pion) photoproduction. 

It is especially interesting to study the scattering matrix S. It is defined on several Riemann sheets of the complex 
energy plane. One finds that poles of S are distributed in a characteristic way on these sheets and can be identified 
with objects like resonances or the cusp. 
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The Lee model was introduced in 1954 by T.D. Lee to study questions concerning renormalization of field theories 
0. Pauli, Glaser and Kallen j^] used this model in order to investigate the appearance of "ghost states", states 
of negative probability density. Ghost states appear when coupling constants are chosen to be large. This lead to the 
discussion of the meaning of an indefinite metric in a Hilbert space. Hohler [ |lO| used the Lee model for a quantum 
mechanical examination of the exponential decay law of unstable particles. The Lee model has not been extensively 



utilized for decades. A good introduction can be found in 11 and |13]. Using the Lee model our main goal was not 



the perfect description of experimental data but to investigate the interplay and meaning of both cusp effect and 
resonance 5*11(1535) at threshold. 

We find that in order to properly treat the threshold effect, rescattering and calculation to all orders of perturbation 
theory is essential. A simple Born approximation at eta threshold is insufficient. We show which Riemann sheets should 
be considered when looking for resonance poles on the complex energy plane. The sheet structure plays an important 
role at threshold. 

From the fact that the experimental data can be roughly described without explicitly using the resonance 5*11(1535), 
we think that the importance of this resonance is commonly overestimated and threshold effects must be taken into 
account. 

In the following we give an introduction to the generalized Lee model. We then discuss the fit results together with 
the experimental data. First we consider the production of eta mesons when pions are scattered by a nucleon. Finally 
we treat eta and pion photoproduction. 

II. A GENERALIZED LEE MODEL 

According to the quantum numbers of the 5ii(1535) resonance, eta production happens mainly in a s-wave channel. 
Since angular momentum and isospin are conserved quantities in strong interaction, we can also restrict our calculation 
to the Hilbert subspace where L = and I = 1/2. Furthermore, we work in the center of momentum system. These 
choices considerably simplify our calculations. 

The particles that appear in our model are the nucleon N, the two resonances Ri and i?2 (namely 5ii(1535) and 
5ii(1650)) and the two mesons Mi and M2 (pion and eta). Taking the quantum numbers in isospin space to be 
/ = 1/2 and = 1/2, the nucleon pion state is actually a superposition of the physical pion states: 

I A^, ^, / = 1/2, h - 1/2) = ^[V2\n)\n+)-\p)\ tt")) . (1) 

We adopt the following notation, so that the isospin does not appear explicitly: 

\N,tt)=\N,ttJ=1/2,I, = 1/2). (2) 

Because we want to make full use of the isospin formalism we have the difficulty of attributing masses to the pion 
and the nucleon since the physical particles {tt°, 7r+) and (n, p) all have different masses. We chose average masses 
= 137 MeV and ttin ^ 939 MeV. 

The mesons and the nucleon are treated as stable particles. However, the resonances 5ii(1535) and 5ii(1650) are 
not stable, because of the decay into nucleon and mesons. This leads to a "dressing" of resonances, e.g. the physical 
resonances are surrounded by a cloud of virtual mesons. In our calculation we use the bare resonance states - not the 
physical ones. Therefore the masses of the resonances, rn_Rj , that appear in the model are bare masses and are used 
as fit parameters in the model. 

The following self adjoint Hamiltonian describes the five particles and their interactions. It consists of two parts, 
the free Hamiltonian and the interaction term: 

H — Ho + Hint , 

so that 

Ho \ Rj) = riiRj I Rj) , 
Ho I N{-k), Mi{k)) = Wi I N{~k),Mi{k)) , (3) 

where Wi is the total energy in the center of momentum (cm.) system. 



Wi = \/mf + fc2 + J m% + fc2. (4) 
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k is the momentum of particles in the cm. system. We distinguish two sorts of interaction among the particles: 
"resonance interaction" and "contact interaction" (see Fig. |l|). In the resonance interaction a meson and a nucleon 
combine to form a resonance, or a resonance disintegrates into a meson and a nucleon 

Rj ^ N,Mi, 



{R, I H,nt I iV(-fc), Ml (fc)) - gf, Fi [Wi) . (5) 

The index I is used for mesons, index j refers to resonances and is the respective coupling constant. FiiWi) is a real 
cut-off form factor which is introduced to assure convergence of expressions that will appear later in the calculations 
for the cross sections. It determines how the interaction between the particles weakens as the kinetic energy rises. 

The contact interaction on the other hand corresponds to a nonresonant interaction, ingoing mesons are scattered 
into outgoing mesons, 

N,Mi ^ N,Mi,, 



{N, Ml (k) I I N, Ml, (fc')> = — Fi, {W{, ) Fi {Wi) . (6) 

The coupling constants gij' for the contact interaction term are symmetric, e.g. gi^f = gv ,i. 

Fig. shows the vertices for these interactions. The cut-off function Fi{Wi) is taken to be the same for the resonance 
and contact interaction terms. This is arbitrary, but it simplifies the calculations. In our computer simulations (see 
sect. we chose Fi{Wi) to be a simple analytical function in order to keep computation time short. In principal, 
Fi can also be deduced from more advanced models by evaluating the matrix elements in Eqs. js] and ^. 

To determine the energy eigenstates of the total Hamiltonian H , we write 

H \N,M,{k))+^W \N,M,{k))+, (7) 

where | N,Mi{k))-\- is a scattering state which satisfies the Lippmann-Schwinger equation for ingoing boundary 
conditions, 

I N, M,{k))+ =1 TV, M,(fc)) - ^— — i^znt I N, M,{k))+ ■ (8) 

Ho — W ~ le 

Mi stands for the incident meson (i ^incident) that at time t = — oo comes in as a plane wave. The Lippmann- 
Schwinger equation can be solved using the following ansatz, 

I iV, M,(fc))+ = E I ^^■) + E / (^') I ^' ^^'(^')) • (9) 

The coefficients a and [3 in Eq. (||) can be obtained by multiplying from the left with {Rj |, {N, Mi \ and using the 
orthonormality of the basis. If we define 



t]^ j d^k Fi{Wi)al{k), (10) 

we get 

^--^^■E^i!^^?- (11) 



where 



am = 6{k - h)6i,, + • J2 ^I'Gvi , (12) 

W — Wi + te ^ 
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Gin can be considered as an effective coupling constant. It is energy dependent and symmetric, Gn = Giii{W) = Gw . 
After inserting ( |T2| ) into ( |lO| ) one obtains: 

Y^hivt], = 6,^iF,{W) (14) 
/' 

with 

hw = - HiGiv (15) 

and 

Eq. (|lj) is a system of linear equations for t\ which can be solved by standard mathematical methods. We note that 
in the derivation of Eq. ( p^ ) the number of different resonances and mesons in the model is arbitrary. If for example 
n different mesons are involved in the scattering process, a n x n matrix needs to be inverted. When solving for ^ 
one encounters an important function which we call hdet- It is the determinant of hw in Eq. ( p^ ) and contains most 
of the information of the scattering process. It is a function of Gw and Hi. As will be shown h^et as well as Hi are 
meromorphic functions of W . They are defined on several Riemann sheets. 
Once the are determined the scattering amplitude Tfi can be calculated. 



Tf, = {N,Mf(kf) I H,^t I iV,Af,(fc,)>+ = -Ff{W)l^Gfv t], , (17) 

v 

where / and i stand for outgoing and incoming mesons respectively. 



A. An example 

In order to familiarize the reader with the Lee model, we present here a simple and instructive example that 
shows the principal features of the general Lec model. We will calculate scattering amplitudes, cross sections and the 
scattering matrix, which will all appear in their familiar textbook forms. 

We consider the case where the two mesons, tt and 77, couple only to one resonance Ri (resonance interaction) 
without contact interactions. The coefficients a and /? of the scattering state are then found to be 

a\{k)^5{k~k:)5u-p^ g'^^^.FiiWi) (18) 

f3' = -9lMW)^, (19) 

n-det 



where /ij„* is the determinant mentioned before: 



hdet = hdet{W) = mR-W- Y^ig^ifHi . (20) 

We immediately obtain expressions for the scattering amplitude Tfi and the scattering cross section 

Tf, = {N, Mfikf) I H,„t I N, M,{ki))+ (21) 
„fl Ff{W)FdW) 

^det 

Using Fermi's rule we obtain the cross section 

da 2tt 1 



-9^ tI^^^J±l^ . (22) 



^ I j^n 



\Tf,\'g{W), (23) 



where jin is the incoming flux and g{W) is the phase space density of the outgoing state. 



4 



After takinga closer look at the Hi we can write the scattering cross section in a very convenient form. We take 
Hi as in Eq. (16) and make use of the identity 



P- 



1 



± iTT 6{Wl - W) , 



Wi-WTie Wi-W 
where P stands for the principal value. Let us write 

Hi = Ai+iTi/2, 

where we set 

§ = (55)2 F?{W) ki"^ e{W -mi- m^) , 



Ai = 4. {gtf P / dWi 

J mi+rriN 



Ff{Wi) , ui E_ 



N 



Wi-W Wi ' 



(24) 
(25) 

(26) 
(27) 



and uji 



kf and E, 



N 



I'M 



+ k?. 



Ti turns out to be the partial width of the resonance for the decay 



channel into meson I and the nucleon. The step function Q in the expression for Ti leads to a zero partial width below 
threshold of meson /, and A; can be interpreted as the mass shift of the resonance due to its interaction with the 
mesons. 

This can be clearly seen when we look at the expression for the scattering cross section that we finally obtain: 



da 



1 



r.r//4 



l-det 



(28) 



r.r//4 



A, 



1^)2 + (r J2 + r„/2)2 



The cross section has the well known Breit-Wigner form of a resonance. We note, that the total Breit-Wigner width 
of the resonance, Ttot, is the sum of the partial widths, F^r + T,,. The physical mass of the resonance, mn — A^r — A^, is 
composed of the bare mass of the resonance ma and the "mass shifts" A^r and A^. The physical mass of the resonance 
is lowered by the coupling to the mesons. The partial "mass shifts" of the resonance. A;, add up to give a total mass 
shift. It is interesting to note that the same function Fi that cuts off the interaction between the particles also is the 
cut-off for the width Fj. Fi is supposed to be a smooth function that gradually falls off as W increases. We then reckon 
that the widths F; grow proportionally to ki{W) at their respective Z-meson threshold when W is increased. This is 
in agreement with general scattering theory that predicts the same behavior for s-wave resonances. 
The S-matrix is a 2 x 2 matrix which corresponds to the two reaction channels. One finds 



S„ = 



iF^r 




n-det 


hdet 


^•\/F7rF^ 


1 

n-det / 


hdet 



So is unitary and the off-diagonal elements are equal which corresponds to time inversion symmetry. From So OIIG CSyll 
easily determine the scattering phases by comparing it to the textbook form of a two channel S-matrix: 



So = 



2i(5i 



„2iS- 



2 ,*(«l+<52) 



where ^ is the inelasticity and S1/2 are the elastic scattering phases. 



B. The Form Factor F, 



The cut-off functions Fi are constructed so that at low energies their values are constant. With increasing energy 
the interaction between the particles gradually decreases. It seems reasonable that at the cm. energy of about two 
nucleon masses the interaction should be strongly suppressed. 
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Although there are many different cut-off functions which fulfiU these features it turns out that in general they 
produce similar results. We selected a cut-off function simple enough so that the integrals Hi could be calculated 
analytically. Explicitly we fixed the cut-off function such that 




where C is a cut-off parameter. C was chosen to be 2000 MeV, so the interaction of the meson with the nucleon stops 
abruptly for Wi > C. We now show that Hi (as a complex function) is defined on two Riemann sheets. 

Let us consider W in Eq. ( p9| ) as a complex variable. Hi is then two sheeted and has a cut from -I- tojv to C. The 
easiest way to see this is by approaching the real axis from above and below and making use of the identity (|2j). The 
continuation of Hi into its second sheet can be done by contour deformation of the integral and applying Cauchy's 
residue theorem. Explicitly we find: 



HiiW) \i.sheet =2V2^iVC-m-mN (30) 

\JC — mi — niN 



— 2\/2mi \/ mi + m^r — W a rrta n ( ■ ) 

^/mi + — W 



Hl{W) \2r^d,heet = HiiyV) Ipts/ieet + 2Try/2^iy^ mi + mN - W . 

Hi diverges aXW — C because of the step of the cut-off at C. In our simulations we always stayed well below C. 
Under the condition C — mi — m^ > > | mi + m^r — | , we can approximate 



Hi{W) \2'-dsheet ~ 2\/2miy/C - mi - mN + TT V2to/ y/ mi + m^ - W . 



The expression ^/mi + m^ — W is directly related to the partial width of the resonance F; as defined in sect. [I A. 
At threshold, e.g. mi + mN = W, y/mi + mN — W becomes imaginary and grows like the meson momentum ki with 
increasing energy W . 

C. The S-matrix 

We can get a better understanding of the physics underlying the scattering and production processes by studying 
the S-matrix So- We will show in the following how poles of the scattering matrix can be identified with resonance 
peaks or cusp peaks in physical cross sections. It will become clear that the distribution of poles and their location on 
different Riemann sheets gives valuable information on the physics involved. We can take full advantage of the fact 
that the Lee model is analytically solvable. Important features of the S-matrix are determined by general physical laws 
like unitarity, causality, time reversal invariance and symmetries of the interactions. A very good introduction to the 
subject can be found in the books by Nussenzveig and Bohm In the following we consider the S-matrix as a 

complex function of the energy W. S is then found to be meromorphic and is defined on several Riemann sheets. It is 
mainly determined by the determinant hdet, which is the same for all channels, e.g. in our case N(tt, 7r)7V, A'^(7r, ri)N, 



N{ri, ri)N and N{ri, 7r)7V. For our example in sect. II A, we have 



hdet = mR~ gl^H, ~ g^.^H^ ~ W , (31) 

i.e. hdet is a four sheeted function, because each Hi is defined on two sheets. There are two cuts located on the real 
energy axis corresponding to the two meson thresholds. They start at the respective meson production thresholds 
which are also the branch points for the Riemann sheets. Fig. ^ shows those four Riemann sheets of hdet, the branch 
points, and the cuts along the real axis. Whenever a cut is crossed one enters another Riemann sheet. Starting from 
the first sheet one moves into the second or third sheet depending on whether the cut was passed below or above 
eta threshold. Passing a sheet twice, brings you back to the original sheet. The fourth sheet can be reached from the 
first sheet by passing the cut twice, once below and once above eta threshold. We find poles in the complex plane 
that belong to the resonance, called resonance poles. The coordinates of the S-matrix poles are functions of the 
coupling constants. In a typical case when the coupling constants are small, the resonance poles are found in the 



second and the third sheet close to the bare mass m^. Following the discussion in sect. II A, the pole coordinates will 
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then be practically determined by the partial widths F/ and the mass shifts A; . If the resonance peak of the scattering 
cross section is located above (below) eta threshold, its shape is determined by the pole in the third (second) sheet. 
As the coupling constants are increased the 'virtual' meson cloud is bound more and more tightly to the physical 
resonance. This lowers the physical mass of the resonance (the resonance peak moves towards the pion threshold) and 
increases its decay width. As a consequence the resonance pole moves away from the physical axis and towards the 
pion threshold (see Fig. ||). In general this behavior is the same for poles in the second or third sheet. However, if the 
coupling constant Qi^^ is held small and constant, and only the coupling of the eta meson gf^,^ is increased, the pole 
in the second sheet moves into the fourth sheet as drawn in Fig. ^. It is this pole in the fourth sheet that gives rise to 
the cusp peak at eta threshold. This will be of importance in the following section. 

We briefly remark that if the coupling of the mesons to the resonance becomes strong enough the physical resonance 
finally becomes a bound state. This bound state is then represented by a pole in the first sheet on the real axis below 
pion threshold. Indeed, with increasing coupling constants the pole of the second sheet moves towards (and finally 
onto) the pion threshold where it can pass across the branch point to the first sheet. This behavior of the pole was 
first seen by Hohler ||T^ (see also in detail and [p^). 



III. SCATTERING PHASE ANALYSIS OF PION NUCLEON SCATTERING 



Because we work in the Hilbert subspace corresponding to quantum numbers / = 1/2, J = 1/2 and L = 0, it is 
convenient to use a partial wave phase analysis of the experimental data to determine the parameters of our model. 
There is a variety of scattering phase analysis which differ partially in their results because the evaluation methods 
and the underlying data are often not identical. Among the best known are those from Karlsruhe Helsinki (KH, 
" |18)), Carnegie-Mellon-Berkeley (CMB, (l|], |§) and from the Virginia Polytechnic Institute (VPI, In our 



work we used mainly the VPI analysis by Arndt and Roper which has the convenience of being supported by the 
SAID interactive dial-in program. 

A typical Argand diagram (Tq) of the 6*11 partial wave for pion nucleon scattering is depicted in Fig. |^. The data are 
taken from a VPI solution (Fall 93) Useful information can be obtained by studying the Argand diagram which 
shows how the complex scattering amplitude Tq for pion nucleon scattering changes as the cm. energy is increased. 
At low scattering energies we get Tq — 0. With increasing energy Tq first moves counterclockwise on the 'unitarity' 
circle which implies purely elastic scattering of the pion. At about the "4 o'clock" position, Tq bends sharply inwards. 
This corresponds to the threshold of eta meson production. It is also at this energy that the resonance S'ii(1535) 
is located. The more Tq moves into the center of the unitary circle, the more the scattering is inelastic, in our case 
due to the production of eta mesons and two-pion production. After a small closed loop To continues in a nearly 
perfect half circle. It is those circular patterns in Argand diagrams that indicate the existence of a resonance p3| , p^ . 
In our case it corresponds to the resonance 5*11(1650). However, there is no circle that corresponds to the resonance 
511(1535). This may imply that either there is no resonance or that its circle is heavily deformed by the opening of 
the eta meson channel. We address this question in the next section when we fit two models, one with and the other 
without resonance 5*11(1535) to the experimental values of Tq. 

Fig.^ shows the calculated elastic (ct^) and inelastic (cr^) partial wave scattering cross sections for the partial wave 
5ii which have been obtained by 

'^s = ^\Tq\^ (32) 
a. = ^(l-|l + 2^^o^. 

The elastic cross section shows two peaks at about 1500 MeV and 1700 MeV. They can be linked to the resonances 
5*11(1535) and 5*11(1650) respectively. The peak at 1500 MeV is located at the eta meson production threshold. Its tip 
is pointed; it has a cusp. The formation of the cusp is a threshold effect and is due to unitarity We show in the 
next section that even the entire corresponding peak at 1500 MeV can be completely explained as a threshold effect. 
The peak can be reproduced in a model without the resonance 5*ii(1535). 

In Fig. H we show the elastic and inelastic cross sections together with the data of Clajus and Nefkens ||2^ for 
the total cross section of eta production. From this comparison we find that the inelasticity of the pion scattering is 
mainly due to the production of the eta meson. For higher energies, however, the multi-pion production is no longer 
negligible and consists essentially of two-pion production in the energy range below 2 GeV. 

Therefore, in order to reproduce the experimental data qualitatively, the Lee model should include two-pion pro- 
duction. However, the inherent Tamm-Dancoff approximation only allows one meson to be present at a time. We 
circumvented this by introducing an additional "meson particle" in our model representing a two-pion system. 
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We then solve the Lee model as before, only now with three different types of mesons instead of two. The third 
meson type stands for the pion pair. The vertices in Fig. ^ show the coupling the two-pion pair to the nucleon and 
the resonances. This introduces three additional coupling constants: gi 2TTi 92 2ttj Q^tt^-k- At the two-pion threshold the 
phase space density for the two-pion particle should not increase proportionally to the momentum k like for the pion 
or eta meson at threshold. Instead for a two-pion system the phase space density grows approximately quadratically 
with energy: 

phase space cx (W27r — S^^r — f^-w)^ ■ (33) 
We can take this into account by choosing H2it to be 

^'-^""); . (34, 

A/.,-(-m„ Mtt {W2^ - W -le) 

The cut-off mass, C2tv, was chosen to be 2000 MeV in our calculation. With our method of representing two pious by 
one effective particle we can not describe exclusive two-pion production, but rather use this representation to explain 
the total reaction cross section. 

By taking into account two-pion production we also change the analytic structure of the S'-matrix. Now there is a 
third threshold at 2m7r + ttin and additional Riemann sheets appear. Moreover with our choice of H2t: in Eq. (jsj), 
H2-K is no longer a two sheeted function of W , but has an infinite number of sheets. However, these sheets do not play 
an important role in our discussion, because we are only interested in the sheet structure in the immediate vicinity 
of the eta threshold. The local sheet str uctu re at eta threshold stays unchanged. Also our discussion of the Riemann 
sheets and the motion of poles in sect. II C is still valid. Fig. |^ shows the new Riemann sheet structure. The sheets 



are numbered such that the sheet structure around the eta threshold looks the same as in Fig. |[ 



IV. DETERMINATION OF COUPLING CONSTANTS AND POLES 



A glance at the partial wave cross section in Fig. ^ shows that at pion threshold (far away from the first resonance) 
the cross section undergoes a rapid change with increasing energy. The elastic pion nucleon scattering cross section 
at pion threshold is even larger than the resonance peaks of the resonances 511(1535) and iS'ii(1650). The great 
magnitude of the elastic scattering cross section is theoretically confirmed by low energy theorems for pion scattering 
which make use of chiral symmetry |^5| . Given the nearly constant form factors at threshold as discussed earlier, 
it is impossible to simulate the rapid fall-off of the cross section from pion threshold to an energy of about 1300 
MeV (see Fig. ||). A better form factor for the pion- nucleon interaction would have to change quite rapidly at pion 
threshold. As a consequence we used our model only at energies above 1350 MeV. 



A. Coupling constants 



The free parameters of the model such as the coupling constants and the resonance masses were determined by 
fitting the Lee model to the experimental data extracted from a partial wave analysis. We fitted two generalized Lee 
models, one including the resonance S'ii(1535), the other without. After fitting we compared the fitting results for 
these two models. From our calculations we observe that the peak at eta threshold is almost entirely due to the cusp 
effect. Even without explicitly including the resonance 511(1535) in the model the peak appears. Moreover we find 
that the experimental data we used in our fits can be qualitatively reproduced. The nonresonant coupling to the eta 
meson gives similar results as obtained with a resonance interaction although the model has four fit parameters less 
than the model with resonance 5ii(1535). Fig. ^ depicts the complex scattering amplitude Tq together with our fitted 
theoretical curves. The data have been taken from the VPI solution (SM90) |^l|. A persistent difference between our 
two models is found at eta threshold (1486 MeV). Im Tq as obtained with the resonance 5ii(1535) is much smoother 
near threshold than without the resonance. Fig. ^ shows the cross section for eta production in pion nucleon scattering. 

In Tab. ||the coupling constants are given for two fits, Lee model with and without the resonance 5ii(1535). Because 
our model is quite different from other standard dynamical models, the coupling constants of the Lcc model can not 
be directly compared to the listed standard coupling constants. 

Also within the Lee model the magnitudes of resonance coupling constants cannot be easily compared to the 
magnitude of contact couplings. This may be understood when looking at their different definitions in sect. In 
order to reproduce the strong cusp effect and the strong production of eta mesons, the overall coupling to the eta 
mesons has to be strong. In one model the eta meson couples strongly to the resonance 5ii(1535) (coupling constant 
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gi^jj) in the other model it couples strongly by contact interaction {gri,-^)- In both models the resonance S'ii(1650) 
couples very weakly to the eta meson while it couples strongly to the pion. 



B. Poles on the Riemann sheets 



One can get more information about the analytic structure of the T-matrix and the peaks of the cross section in 
Fig. H by looking at the pole distribution on the Riemann sheets. Tab. |l| gives the coordinates of the poles for the 
two different models and Fig. |l^ shows their locations in the complex plane. We note that the resonance S'ii(1650) 
in both models has poles on the 2'"^ and the 3'"'' sheets, as is to be expected from sect. 11 C . The coordinates of the 



poles agree well with the mass and the width of the resonance 5*11(1650). In both models we find a pole on the 4*'* 
sheet close to the eta threshold which is responsible for the cusp peak. This pole has not previously been discussed. 
No additional resonance pole in the 3'''^ or 2"'' sheet is needed to account for the cusp peak. We can trace back the 



origin of this pole by numerically lowering the respective eta meson coupling constants as described in sect. |IC|. The 
pole then moves back to its original starting point, which is different for the two models. For the case of the model 
with resonance S'ii(1535) the pole is originally a resonance pole in the 2"^^ sheet that has been pushed into the 4*'' 
sheet by the strong coupling to the eta meson, whereas in the model without resonance 5*11(1535) the pole always 
stays in the 4*'' sheet. 

Another pole in the 3'''' sheet for the case of the model with resonance 5*11(1535) is the resonance pole of the 
resonance 5ii(1535). 

Finally we have compared our direct determination of the poles with the speed plot technique suggested by Hohler 
0. The pole of the resonance 5ii(1650) in the 3'"'' sheet is very well reproduced by the speed technique and found 
as (1666-811) MeV in our model with 5*11(1535) and as (1671-571) MeV in our model without 5ii(1535). In addition 
in both our models we see a very narrow spike with a width of a few MeV at the 77 threshold which corresponds to 
the cusp. The resonance pole of the 5*ii(1535) in our model does not show up using the speed technique. It is entirely 
covered by the cusp effect. 

V. ETA AND PION PHOTOPRODUCTION 

In order to check the consistency of our model and to further study the role of the resonance 5*11(1535), we studied 
eta and pion photoproduction. Since the photon coupling is weak, we use perturbation theory for the photoproduction. 
The strong interaction is still fully accounted for as we can employ the full analytical solution of the strong interacting 
particles discussed in the preceding sections. Fig. ^ shows how the photon couples to the hadrons via resonance and 
contact interactions. The matrix elements are defined as follows: 

{R, I H, I Ni-i^),j{,7),e) ^^d-e, (35) 



(iV(-i7),7(i;),e I I N{-k),Miik)) - FiiWi) a ■ e, (36) 



where j is again the index for the resonances, and gi^^ are the coupling constants for resonance and contact terms, 
respectively. The momentum and the polarization vector of the photon are denoted v and e, and a is the spin of the 
nucleon. 

To first order approximation of perturbation theory, the scattering matrix Tfi is then 

Tf, = 4N, Mf{k) \H^\N, 7(i7), e) , (37) 
and -{N, Mf \ is the outgoing scattering state as calculated in the previous section, 

_(iV, Mf 1= Pi (R, I + E / (^) {N{~k),Mi{k) I . (38) 

The coefficients a and /3 are given in Eqs. (^ and (|l|). As before, we consider the two models with and without 
the resonance 5ii(1535), with their respective parameters and coupling constants given in Tab. | for pion scattering. 



We fitted to the eta photoproduction cross section and pion photoproduction amplitude i?o+ (see Tab. III). The 
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experimental data are taken from Wilhclm and Krusche for eta photoproduction and for pion photoproduction 
we used the VPI solution (Sp 95) |2l[]. The results are shown in Figs. |lj and |l^. The fit without resonance is not as 
good as the fit with the resonance. This is partially due to the smaller number of fitting parameters (five parameters 
less). However, neither model can well describe the dip of Im(£^o+) at 1600 MeV. A more realistic description of the 
background would be necessary to better describe the data below rj threshold. This would also improve the result for 
eta photoproduction for the case without the 5ii(1535) resonance. 

VI. CONCLUSION 

We have constructed a generalized Lee model for pion scattering and eta production in the s-wave channel. Using 
only the Hilbert subspace corresponding to the quantum numbers of the partial wave 5*11, the particles that appear in 
the model are the resonances 5*11(1535) and S'ii(1650), the nucleon, the pion and the eta meson. These five particles 
interact via a resonant and a nonresonant interaction. The absence of antiparticles makes it a nonrelativistic model. 
However, it has the advantage of being fully unitary and analytically solvable. Eta and pion photoproduction were 
additionally calculated using a perturbation approach. The coupling constants and other parameters of the model 
were determined by fitting the model to partial wave analyses and cross sections. The experimental data we used 
were total cross sections for eta meson production in pion scattering and eta photoproduction [p6| . In addition we 
employed a partial wave analysis of pion nucleon scattering and £^o+ amplitudes for pion photoproduction Q , |^ . 

We were interested in studying the cusp effect and in obtaining a better understanding of the interplay of the 
resonance 511(1535) and the cusp at eta threshold. This is particularly of importance when determining why the 
neighboring resonances 5ii(1535) and 5ii(1650), having the same quantum numbers, show quite different behavior: 
5ii(1535) couples strongly to the eta meson while the resonance 5ii(1650) couples to it only very weakly. Following 
Hohler's statement Q that partial wave analyses show no evident signature for the resonance 5ii(1535), we investigated 
a Lee model with and without the resonance 5ii(1535). 

We found that using the Lee model the experimental scattering data for various scattering channels could be qual- 
itatively reproduced without introducing the 5ii(1535) resonance. Therefore we think that the resonance S'ii(1535) 
is weaker and less important than generally accepted. 

The S-matrix, together with its Riemann sheet structure, was studied. Eta threshold is the branch point where the 
four Riemann sheets meet that were of importance in our discussion. We found a specific distribution of poles on these 
sheets that determines the scattering amplitude on the physical axis. Poles could be attributed to resonance peak or 
the cusp peak. More precisely a pole in the 4*'' sheet that until now has never been taken into account, determines 
the shape of the cusp. 

When doing calculations close to eta threshold, rescattering and full calculation to all orders are important. 
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TABLE I. Coupling constants and masses of our two models. 











Model with resonance 


511(1535) 










couplings 




91,. 


92,r, 


52. ir 5l,27r 


92,2tt 


9v,v 


9ri-'^ 


QtT ,7V 


Qtt ,'2-K 


values 


0.24 


0.12 


0.01 


0.25 0.0023 


0.062 


2.3x10"* 


0.01 


0.0031 


0.024 



mm = 1616 MeV, mR2 = 1713 MeV. 











Model without resonance S'ii(1535) 



5^ 
-0.24 



5l,27r 




92,2. 

0.052 



couplings 
values 



5»;,i7 
-0.052 



5'7.T 

0.045 



p7r,7r 

0.0083 



57r,27r 

-0.0098 



mR2 = 1704 MeV. 



TABLE II. Pole positions on the Riemann sheets for our two models. 



model 


5ii(1650) 




Cusp/5ii(1535) 




2. sheet 


3. sheet 


3. sheet 4. sheet 


with 511(1535) 


(1655;-111) 


(1662;-87) 


(1501,-61) (1528;35) 


no 5*11(1535) 


(1652;-90) 


(1670;-60) 


(1458;57) 



TABLE III. Electromagnetic coupling constants for eta and pion photoproduction. 



coupling constants 


5^7 


52,7 


9t^,^ 


517.7 


with 5ii(1535) 


-2.18 


1.58 


0.0027 


-0.0011 


without 5ii(1535) 





1.3 


0.0024 


-0.0058 
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FIG. 1. The vertices for the resonance and the contact interaction. The indices I and I' stand for the meson type, the index 
j refers to the resonances. 
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FIG. 2. The four Riemann sheets for the pion and eta meson and the paths from one sheet to the next. 
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FIG. 3. Pole movement in the 3'''' and the 2""* (4"^) sheet (see text). 
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Argand Diagram 
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FIG. 4. Argand diagram for the ttN partial wave Sn. Data from VPI [ pl[ . 
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Partial Wave Cross Sections 
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FIG. 5. Scattering cross sections for the partial wave Sn, calculated from the partial wave analysis by Arndt and Roper 

( ) elastic cross section (Js for nN — > ttN, ( ) inelastic cross section a"r. The data points are the cross sections for 

eta production p3|. 
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FIG. 6. Vertices for the two-pion system. 
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FIG. 7. Riemann sheets when two-pion production is taken into account. 
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FIG. 9. The cross section for eta production in pion nucleon scattering. The solid and dashed lines show the model 
calculations with and without resonance 5*11(1535). Experimental data are taken from Nefkens |23[. 
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FIG. 10. The pole positions on the Riemann sheets for our two models. (*) model with 5'ii(1535), (□) model without 
511(1535). 
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FIG. 12. The total eta photoproduction cross section with resonance S'ii(1535) ( — ) and without 5ii(1535) ( ). 
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